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Abstract 

We consider a scenario (supported by some lattice results) in which a [/(l)-breaking 
condensate survives across the chiral transition in QCD. This scenario has important 
consequences for the pseudoscalar-meson sector, which can be studied using an effective 
Lagrangian model. In particular, generalizing the results obtained in two previous papers, 
where the effects on the radiative decays r], rf — > 77 were studied, in this paper we study 
the effects of the U(l) chiral condensate on the strong decays of the "light" pseudoscalar 
mesons, i.e., 77,7/ — > 3tt°; T],T]' — > 7r + 7r~7r°; 7/ — > r/7r°7r°; 7/ — > r]7T + 7T~; and also on the 
strong decays of an exotic ("heavy") S*?7(3)-singlet pseudoscalar state rjx, predicted by 
the model. 
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1. Introduction 



It is well known that the QCD vacuum has a very complicated structure, characterized 
by some non-trivial local (or also non-local) condensates, whose behaviour as a function 
of the temperature T also characterizes the phase structure of the theory. 
For example, a phase transition which occurs in QCD at a finite temperature T c ^ is the 
restoration of the SU(L) Cg> SU (L) chiral symmetry (in association with L = 2, 3 massless 
quarks), which for T < T c h is broken spontaneously by the non-zero value of the so-called 
11 chiral condensate", i.e., (qq) = X^=i(<Mi) PQ- But QCD with L massless quarks has also 
(at least at the classical level) a (7(1) axial symmetry [21 E]. This symmetry is broken 
by an anomaly at the quantum level, which in the " Witten-Veneziano mechanism" [HE] 
plays a fundamental role (via the so-called "topological susceptibility") in explaining the 
large mass of the if meson. The role of the U(l) axial symmetry for the finite temperature 
phase structure has been so far not well clarified. One expects that, above a certain critical 
temperature also the £7(1) axial symmetry will be (effectively) restored but it is 

still unclear whether Tf/(i) has or has not something to do with T c h. 

In this paper we re-consider a scenario (which was originally proposed in Refs. [61 
[71 El [9] and elaborated in Refs. [TUl CHI H2] , and which seems to be supported by some 
lattice results on the so-called "chiral susceptibilities" [131 US]) in which a new £7(1)- 
breaking condensate survives across the chiral transition at T c h, staying different from zero 
up to a temperature Tun.) > T ch . Tun.) is, therefore, the temperature at which the £7(1) 
axial symmetry is (effectively) restored, meaning that, for T > Tun), there are no £7(1)- 
breaking condensates. The new £7(1) chiral condensate has the form Cu{\) = (Ou(i)), 
where, for a theory with L light quark flavours, Oum is a 2L-fermion local operator that 
has the chiral transformation properties of [31 [161 D2J 

Outi) ~ det(q sR q tL ) + det(q sL q tR ), (1.1) 

st st 

where s, t — 1, . . . , L are flavour indices. The colour indices [not explicitly indicated in Eq. 
f ll.il) ] are arranged in such a way that: i) Oum is a colour singlet, and ii) Cjjm = (Ou(i)) 
is a genuine 2L-fermion condensate, i.e., it has no disconnected part proportional to some 
power of the quark-antiquark chiral condensate (qq): the explicit form of the condensate 



Throughout this paper we use the following notations for the left-handed and right-handed quark 



fields: q L . R = ±(1 ± 75)9, with 75 = -i7°7 1 7 2 7 3 
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for the cases L = 2 and L = 3 is discussed in detail in Appendix A (see also Refs. 

This scenario has important consequences for the pseudoscalar-meson sector. The 
low-energy dynamics of the pseudoscalar mesons, including the effects due to the anomaly, 
the qq chiral condensate and the new U(l) chiral condensate, can be described, in the 
limit of large number N of colours, and expanding to the first order in the light quark 
masses, by an effective Lagrangian written in terms of the topological charge density 
Q, the mesonic field Uij ~ qjRqn (up to a multiplicative constant) and the new field 
variable A ~ det (q S RqtL) (up to a multiplicative constant), associated with the new U(l) 
condensate El [10] : 

C(U, U\ X, X\ Q) = i Tr^l/Wrf) + -d^Xd^X^ 
-V(U, U\ X, A f ) + 1 -oj x Q Tr(ln U - In U ] ) 
+1(1 - Ux)Q(hxX - In At) + (1.2) 
where the potential term V(U, U\X,X<) has the form: 

V(U, U\ A, At) = £ Tt[(rfu - p.I) 2 ] + tf(X*X - Px f 

Bm : Tr(MC7 + Aftf/t) _ J*[det(U)X* + det(U^)X]. (1.3) 



2^2 2^2 

M = diag(mi, . . . , mi) is the quark mass matrix and A is the topological susceptibility 
in the pure Yang-Mills (YM) theory. (This Lagrangian generalizes the one originally 
proposed in Refs. [HJ EEJ EUJ ED |22] , which included only the effects due to the anomaly 
and the qq chiral condensate.) All the parameters appearing in the Lagrangian must 
be considered as functions of the physical temperature T. In particular, the parameters 
p n and px determine the expectation values (U) and (A) and so they are responsible 
respectively for the behaviour of the theory across the SU(L) <g> SU(L) and the U(l) 
chiral phase transitions, as follows: 

1 2 

Pn\T<T ch = 2 F - >0 ' P*\T>T ek <0, 

Px\t<t u(1) = ~F X > 0, Px\t>t u(1) < 0. (1.4) 

The parameter F n is the well-known pion decay constant, while the parameter Fx is 
related to the new U(l) axial condensate. Indeed, from Eq. fll.4p . px = \F X > 
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for T < 7j7(i), and therefore, from Eq. fll.3j) . (X) = F x /\/2 ^ 0. Remembering that 
X ~ det (q S RqtL), U P to a multiplicative constant, we find that Fx is proportional to 
the new 2L-fermion condensate = (Oum) introduced above. In the same way, the 

pion decay constant F n , which controls the breaking of the SU(L) £g> SU(L) symmetry, is 
related to the qq chiral condensate by a simple and well-known proportionality relation 
(see Refs. |H ID] and references therein): (qiqi)T<T ch — —\B m F^. (Moreover, in the 
simple case of L light quarks with the same mass m, m 2 NS = mB m /F n is the squared 
mass of the non-singlet pseudoscalar mesons and one gets the well-known Gell-Mann- 
Oakes-Renner relation: m 2 NS F 2 ~ —2m(qiqi) T<Tch .) It is not possible to find, in a simple 
way, the analogous relation between F x and the new condensate Cf/(i) = 

However, as was shown in two previous papers [TU information on the quantity 
Fx (i.e., on the new U(l) chiral condensate, to which it is related) can be derived, in 
the realistic case of L = 3 light quarks with non-zero masses m u , rrid and m s , from the 
study of the radiative decays of the pseudoscalar mesons r] and rj into two photons. A 
first comparison of the results with the experimental data has been performed and it is 
encouraging, pointing towards some evidence for a non-zero U(l) axial condensate. 
The following decay rates are derived [TT], [T2] : 

r(??^77) = 192^2 (cosy5+— ^siny ) . (1.5) 



„2™3 
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r(n' — >■ 77) = ^-tt cos lo — sin oj 1 , (1.6) 

v/ " ; 192vr 3 F2 I F v > I K ' 

where a = e 2 /47r ~ 1/137 is the fine-structure constant. Here F v > is defined as follows: 



Ft = y/F* + 3F%, (1.7) 

and can be identified with the 7/ decay constant in the chiral limit of zero quark masses. 
Moreover, (p is a mixing angle, which can be related to the masses of the quarks m u , m^, 
m s , and therefore to the masses of the octet mesons, by the following relation: 

y/2 _ F F I 

where: ml = 2Bm and m 2 = |-B(m + 2m s ), with: B = ^f - and m = m "+ md . 

If one puts F x = 0, i.e., if one neglects the new U(l) chiral condensate, the expressions 
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written above reduce to the corresponding ones derived in Ref. [23] using an effective 
Lagrangian which includes only the usual qq chiral condensate. Using the experimental 
values for the various quantities which appear in Eqs. ( 11.51) and 



61), one can extract 



the following values for the quantity Fx and for the mixing angle (p\l 

F x = 24(7) MeV, <p = 17(2)°, (1.9) 

and these values are perfectly consistent with the relation (11.81) for the mixing angle, if 
one uses for the pure-YM topological susceptibility the estimate A = (180 ± 5 MeV) 4 , 
obtained from lattice simulations 



In the section 3 of this paper, continuing the work started in Refs. [Tlj [12], we 
shall study the effects of the U(l) chiral condensate on the strong decays of the "light" 
pseudoscalar mesons, i.e., 77,77' — > 3tt°; T],T]' — > 7r + 7r~7r°; 7/ — > f]TT 7T°; 77' — > 777r + 7r~; 
and also on the strong decays of an exotic ("heavy") S'L r (3)-singlet pseudoscalar state 
r/x, predicted by the model: r/x — > 37r°; r/x —> n + n^n ; r/x — > 777r°7r°; i]x —> r]7Y + n^; 
i]x — > r]'7r 7t ; t]x — > r]'7i + 7T~] r/x — > 3?7, 37^', r)7]r]' , r]7]'i]'. In particular, in the case of the 
exotic particle t]x, we shall find some relations between its mass and its decay widths, 
which in principle might be useful to identify a possible candidate for this particle. 

For the benefit of the reader, we shall start, in section 2, by resuming the main 
results, obtained in the original papers [6j [HI [10], concerning the mass spectrum of the 
Chiral Effective Lagrangian (ll.2p -f fl~3|) . for temperatures T < T^: in this paper we shall 
consider the case T = only. 



indeed, the original values reported in Refs. [lU H2] were: F x = 27(9) MeV and (p = 16(3)°. The 
values reported in Eq. (|1.9|) (which are, anyhow, consistent with the original values within the errors) 
have been obtained using the updated experimental values of the Particle Data Group 24J (in particular: 
r cxp (?7 77) = 0.51(3) keV and r exp (?/ -> 77) = 4.31(36) keV; moreover we use: = 92.2(4) MeV, 
~ 134.98 MeV, m, q ~ 547.85 MeV, m v > ^ 957.78 MeV). 
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2. Mass spectrum and new parameters of the Chiral 
Effective Lagrangian 

Let us consider the Lagrangian ( 11.21) . where the field variable Q(x) has been integrated 
out: 

c(u, u\ x, xt) = I Ti(d^ud^) + i^iyi 1 - V(U, U\ X, xt) 

+-A [ Wl TT(lnU -InU 1 ) + (1 - Wi)(lnX - lnX f )] 2 . (2.1) 



2.1. Mass spectrum at T = for a generic L (in the chiral limit) 

At T = both SU(L) eg) SU(L) and U(1)a symmetries are broken. 

Following Ref. [18], we can eliminate the redundant (having much larger masses) scalar 
fields of the linear a-type model by taking the limit A^ — > oo and A^ — > oo. In this limit 
the potential term gives the following constraints: 

^U=\fI-\- X*X = \f%. (2.2) 

We are thus left with a non-linear chiral effective model, in which the field U has the 
form: 

L 2 -l 

U=\l exp { ^$ }, $ = V n a r a + ^=1, (2.3) 



where r a (a = 1, • • • , L 2 — 1) are the generators of SU(L) (Tr(r ) = 0) in the fundamental 
representation, with normalization Tr(r a T5) = 5 a t,, and 7r a (a = 1, • ■ • , L 2 — 1) are the 
non-singlet meson fields, while S n is the usual quark-antiquark SU (L) singlet field: 

L 

Sir~iy*] (qiLqm - (2.4) 
1=1 

And similarly the field X has the form: 

X=\F~Fxew{ l 4^Sx) (2.0) 



2 F 
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where Sx is an exotic singlet field, with the following quark content: 



S x ~ i[det(q sL q t R) - det(q sR q tL )}- 

st st 



(2.6) 



Substituting Eqs. (12. 3 P and (12. 5 P into Eq. (12. ip and taking only the quadratic part of the 
Lagrangian, we obtain: 



~^7r ^7r a + -d.S^Sn + \d,S x d»S x - \ (E^ t ^J n a 



where: 



1 (y/2L V2 q \ 1 
2 C [^-^ ] ~- A 



<2L V2 
-p— + — (1 - ujijbx 



-i 2 



2 B m 



In the chiral limit, sup rrii — > 0, Eq. (12. 7p reduces to 



-mi. 



£ 2 



2 



/2L >/2 \ 1 
' "FT " ^ " 2^ 



(2.7) 



(2.8) 



(2.9) 



In this case the L 2 — 1 non-singlet fields are massless: they are the Goldstone bosons 
coming from the breaking of the SU(L) £g> SU(L) symmetry down to SU(L)y. Instead, 
the two singlet fields S n and Sx are mixed with the following squared mass matrix: 



2L(Auii+c) 2 v / L[Aji(1-ui)-c] ' 

Fl Fit Fx 

2\/L{Auj 1 (1-uj 1 )-c] 2[A(1-lj 1 ) 2 +c] 
F^F X 



(2.10) 



r x 



The eigenvalues of this matrix are: 



where: 

2A[F?(1 - uj,Y + LF x uf] + 2c(f? + LF\) 



Z 



L 



7T A 



4LAc 



Ql = T^T- (2-12) 

Making use of the following iV-dependences of the relevant quantities in the limit of large 
number of colours N (see Ref. [6]): 

F 7T = 0{N 1 / 2 )- F x = 0(N 1 ' 2 )- A = 0(1); c = 0{N), (2.13) 

we derive, at the first order in the 1/N expansion (and assuming that c\ ^ 0: see the 
discussion in Appendix B), the following expressions for the two eigenvectors: 



Si = - 7 =L==(f 7T s 7T + Vlf x s x) 

V F n + LF X 



S 2 = 1 (y/LF x S„-F„S x ), (2.14) 

V^TT + Lt X 

with the corresponding eigenvalues: 



The two fields Si and S 2 have the same quantum numbers, but different quark contents: 
the first one (assuming that F n ^> Fx) is prevalently a quark-antiquark singlet S w , 
while the second one is prevalently an exotic 2L-fermion singlet Sx ~ ^[det(g s£ g 4 /j) — 
det(g s ^g t i)]. Both fields are massive in the chiral limit. If we let F x — > in the above- 
reported formulae (i.e., if we neglect the new U(l) axial condensate), then Si —> S n and 
nig — >■ 2LA/F 2 , which is the usual Witten-Veneziano formula for the rf mass in the 
chiral limit |H [5]. On the other side, m| 2 ~ 2c/ F x — > oo for F x — > 0, being c = 0(F X ) 
[Eq. f )2.8p ]. and therefore, in this limit, the field S 2 — > —Sx is "constrained" to be zeroo 



*More rigorously, before taking the limit Fx — > (i.e., X — > 0), one should first take the limit ui\ — > 1, 
so that no singular behaviour arises from the anomalous term in Eqs. (|1.2|) and (|2.1j) and the Lagrangian 
simply reduces, for X — > 0, to the usual Lagrangian of Witten, Di Vecchia, Veneziano et al. It is easy 
to check that, by putting uj\ — 1 in Eqs. (|2.11|) - (|2.12j) and then letting Fx ~ > 0, one recovers the 
same results that one also obtains by simply letting Fx — > in Eqs. (|2.15l) . i.e., m% 1 — > 2LA/F% and 
m| ~ 2c/F x -> oo. 
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In the more general case Fx 7^ 0, which we are considering in this paper, there is a field 
(Si) with a squared mass which vanishes as 0(1/ N) in the large-iV expansion; on the 
contrary, the field S2 has a large mass of order 0(1) in the large- N limit. It is quite easy 
to convince oneself that the particle associated with the field Si is nothing but the particle 
if, which is required by the well-known Witten-Veneziano mechanism for the solution of 
the U(l) problem (see Refs. [HI ED])- In fact, the expression for the U(l) axial current: 

42 = ilTrirfdpU-Ud^+LitfdftX-XdpX*)) 

= -VUd^F^ + VLFxSx), (2.16) 

can be re-written, using the first Eq. ( I2.14p . as: 

J 5 ( 2 = -V2LF Sl d,Si (2.17) 

where : 

F Sl = JfJ+LFx (2.18) 
is nothing but the decay constant of the singlet meson Si, defined as: 

<0|4>)|Si(Pi)> = »V2lFaJV (2.19) 



We remind that, according to the Witten-Veneziano mechanism for the solution of the 
U(l) problem, the r/ mass must satisfy the following relation, known as the Witten- 
Veneziano formula: 

m% = — r . (2.20) 

Using the first Eq. (12. 15p , together with Eq. ( 12.18P , one immediately verifies that the sin- 
glet meson associated with the field Si indeed verifies this relation, i.e., m% = 2LA/Fg . 
For this reason, from now on, the field/particle Si will be denoted as rj', with: 



F v , = F Sl = ^F2 + LF x . (2.21) 

Instead, from now on, we shall use the name rjx to denote the other exotic singlet 
field/particle £2. 
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2.2. Mass spectrum at T = for the realistic L = 3 case 



Let us consider more carefully the realistic case [8], in which there are L = 3 light quark 
flavours, named u, d and s, with masses m u = (1.74-3.3) MeV, = (4.14-5.8) MeV and 
m s = (80 4- 130) MeV [24J, which are small compared to the QCD mass-scale Aqcd ~ 
0.5 GeV. In this case Eq. (I2.3P becomes: 

8 

<& = $> a r a + -^I, (2.22) 

a=l * 

where n a (a = 1, • • ■ , 8) are the pseudoscalar mesons (J p = 0~) of the octet, while S n 
is the quark-antiquark S'6 r (3)-singlet field. Proceeding as in the previous section, but 
making also an expansion up to the first order in the quark masses, we immediately find 
that the fields tti, 7r 2 , 7r 4 , 7r 5 , 7r 6 , tt 7 are already diagonal, with masses: 

< 2 = ml ± =B(m u + m d ), 




<, s = m> K± = B{m u + m s 

4 6 ,7 



m ^67 = m K° Ro = B{m d + m s ), (2.23) 



where B = 

On the contrary the fields 7r 3 , tc 8 , S n , Sx mix together, with the following squared mass 
matrix: 



( 2Bm 



7a SA 
%BA 



\ 







\B{m + 2m s ) 
^B(m-m s ) 




IBA 

^B(rh-m s ) 

TV 

2y / 3[A(l-q;i)m-c] 
FkF x 



\ 



2y / 3[A(l-q;i)hJi-c] 
FkF x 
2{A(l-u n ) 2 +c] 

n / 



(2.24) 



where rh = m "^" md , = |5(2m + m s ) and A = m u — m d . This last parameter A 
measures isospin violations, i.e., the explicit breaking of the SU(2)y symmetry. If we 
neglect the experimentally small violations of the SU(2)y isospin symmetry, i.e., if we 
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put A = in Eq. f)2.24|l FI. the squared mass matrix ( I2.24p simplifies to: 



/C 



(2Bm 



v 





|.B(m + 2m s ) 
^0-B(fh — m s ) 







^B(m - m s ) 

60M+c) , 2 
2y / 3[i(l-u;i)u;i-c] 



\ 


2y / 3[A(l-g;i)h;i-c] 

2[A(l-u, 1 ) 2 +c] 




(2.25) 



/ 



Therefore, in this limit, 7r 3 also becomes diagonal and can be identified with the physical 
state 7T°, with squared mass: 



2Bm = B{m u + ma) = rn\ 



(2.26) 



The fields (tt 3 ,tt 8 , S n , Sx) can be written in terms of the eigenstates (Tr°,T],r)',r)x) as 
follows: 



(2.27) 







/tt°\ 




= C 


V 
rf 


\Sx) 




W 



where Co is the following orthogonal matrix (TTJ [12] : 



/I 







cos <f 





— sin (p 



jrj Sin (f Jr^ COS 



yi ^sin^ 



F fcos</? 



\ 



F V> 

-Ftt 



(2.28) 



As we have already said above, ify = ^ F% + 3F X can be identified with the 7/ decay 
constant in the chiral limit of zero quark masses [Til [12]. Moreover, (p is a mixing angle, 
which can be related to the masses of the quarks m u , m^, m s , and therefore to the masses 
of the octet mesons, by the relation ( 11. 8ft [TT1 [12] . 

The matrix Co has been derived by diagonalizing the squared mass matrix (I2.25P at 
the first order in the quark masses and in 1/N, so neglecting terms behaving as 1/iV 2 , 
m 2 or m/N (and assuming, again, that c\ ^ 0: see the discussion in Appendix B). 



tin the next section, instead, we shall take into account also the small violations of the SU{2)v isospin 
symmetry, by taking A^O. 



11 



Following Refs. [SJ [TBI [23], we have considered the limit in which raj Aq C d <C 1/iV <C 1: 
this particular choice is justified by the fact that the mixing angle, which is of order 
0(mN/ Aqcd), is experimentally small*. The other eigenvalues of the squared mass matrix 
( 12.251) can be easily derived at the first order in the quark masses and in 1/N (in the sense 
explained above): 

(2.29) 
(2.30) 

m 2 . (2.31) 

The physical interpretation of these three states is clear. The state 77 is the eighth pseudo- 
Goldstone bosons of the octet: its mass vanishes with the light quark masses. On the 
contrary, the states rf and r\x have masses which do not vanish with the light quark masses. 
In particular, the state rf has a topological (non-chiral) squared mass term 6A/F 2 ,, which 



m\ 


2 

= -B(m + 2m s ), 




m\, 


6^4 Fl 2 

rf r rf 






2cF* 2A[F> 1 -1)+3F> 1 ] 2 


3F X 




F 2 



vanishes as 1/N in the large- N limit. The state r/x, instead, should be heavier, having a 
normal (non chiral) mesonic mass term_j of order 0(1) in the large- N limit. 
From Eqs. ( 12.231) . ( 12.261) and (12. 291) one immediately derives the well-known G ell-Mann- 
Okubo formula [281 12H] for the squared masses of the octet mesons: 

3mJ + ml = 4m|, (2.32) 

where: m\ = |(m|. ± + m? KQ R0 ) = B(m + m s ). In fact, it is natural to expect that the 
introduction of a new chiral order parameter, which only breaks the U(l) axial symmetry, 
should not modify the mass relations for the octet mesons, such as Eq. (1 2 . 3 2 j) . which only 
derive from the breaking of SU(3) <S> SU(3) down to SU(3)v- 

Considering also the squared mass (12.301) of the rj', one immediately derives the following 
interesting relation (with m\ defined as in Eq. (I2.32p ) [S]: 

1 + 3^ m\, + m 2 ri - 2m\ = ^. (2.33) 



tin the literature, also other possibilities have been studied. For example, Leutwyler in Ref. [26] 
considers iti/Aqcd and 1/N to be of the same order, and Witten in Ref. [TH] studies also the opposite 
case, i.e., tuN/Kqcd ^> 1. 

§See Ref. [37] for a detailed discussion of hadrons and their masses in the framework of the 1/N 
expansion. 
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This is nothing but a generalization of the usual Witten-Veneziano formula for the rf mass 
(including non-zero quark masses), with a correction which only depends on the parameter 
Fx (which, as we have already said in the Introduction, is essentially proportional to the 
new U(l) axial condensate), but not on the other unknown parameters of the model (u\, 
Ci). From Eq. (I2.33p . using the known values for the meson masses, the pion decay 
constant F n and the pure-gauge topological susceptibility A, one can derive the following 
upper limit for the parameter Fx- \Fx\ < 20 MeV [EJ [TO] . 

Finally, we can derive an anologous relation involving also the squared mass of the exotic 
state rjx- By taking the trace of the squared mass matrix (12.24ft . using the relations (I2.23jl . 
together with m = m " ^ d , m\ = ~B(2m + m s ) and m 2 K = B{m + m s ), one obtains: 

Tr[/C] = mlo + m 2 + m^, + m 2 x 

„_ 2„, , 2 6(Au 2 + c) 2A(\-uj 1 ) 2 + 2c 

6 TV ?X 



, 6(Aco 2 + c) 2A(1 -u 1 ) 2 + 2c 

mlo + 2m\ + 1 ; 2 ' + — * =^ , (2.34) 

r„ r x 



from which, re-ordering, one finally gets: 

2 2 2 2 2cF v> 2A\F 2 (l - U!) 2 + 3F 2 u 2 ] 
<c + <• +m 2 -2m 2 K = -^ 2 -+ [ ^ -JL ^1. (2.35) 

Unfortunately, this expression depends upon all the unknown parameters of the model 
(Fx, U)i, Ci) and, therefore, we cannot use it to obtain a direct estimate of the mass of the 
particle rjx- However, in the next section we shall find some relations between its mass 
and its decay widths, which in principle might be useful to identify a possible candidate 
for this particle. 
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3. The strong decays of the pseudoscalar mesons 

71,*/, VX 



In this section we shall study the strong decays of pseudoscalar mesons, using the Chiral 
Effective Lagrangian which we have discussed above. 

First we observe that the strong decays of a pseudoscalar meson into two pseudoscalar 
mesons are trivially forbidden by parity conservation. In fact, in terms of the Chiral 
Effective Lagrangian (12. ip . one easily verifies that it is invariant under the following field 
transformation: 

U^U\ X^X\ Q->-Q, (3.1) 

which is nothing but the parity transformation for the fields [provided one also transforms 
the space-time coordinates as x = (x°, x) — > xp — (x°, —x)]. In terms of the meson fields 
7r a , S n , Sx, defined in Eqs. (12. 3p and (12. 5p . it corresponds to: 

n a -> -n a , S n -> -S w , S x -> -S x . (3.2) 

Therefore, terms with an odd number of meson fields (which are not parity invariant) 
necessarily vanish. In particular, operators with three pseudoscalar meson fields are absent 
and therefore the strong decays of a pseudoscalar meson into two pseudoscalar mesons 
are forbidden. 

On the contrary, the strong decays of a pseudoscalar meson into three pseudoscalar mesons, 
being induced by parity-invariant four-meson operators, are allowed and we shall devote 
the rest of this section to a detailed discussions of these decays. 



3.1. The four— meson Lagrangian 



In order to study the strong decays of 77, rj', r] X into three pseudoscalar mesons, we have 
to isolate the four-meson operators in the Lagrangian (12.11) . when expanding the fields 
(12.31) and (12.51) in powers of the meson fields. We thus obtain the following four-meson 
Lagrangian: 



+ 



4F 2 
B 



Tr 



3 



4F1 



d tl S x d»S x + ~S x ns x 



Tr [M$ 4 ] + 



Vs 1 
^ " T~ x Sx 



(3.3) 
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where, as usual: B = ^s 2 -, c 



2Kr' 



J* (Ex \ (fA 3 

By making an integration by parts and using the usual identitites for the SU (3) generators, 
we can re-write the first term in the r.h.s. of Eq. (13.31) as (apart from total derivatives): 



5C 



4i* 
1 



Tr 



Tr 



4 

3°" 



(3.4) 



where / a ft c are the structure constants of SU(3), defined as: [r a ,Tb] = i\f2f abc T c , with 
Tr(r a r fe ) = <5 ab . It is easy to see that this term gives contributions only to decays into 
charged pions, whose fields are 7r = 7ri ^ 7r2 . 

Concerning the second term in the r.h.s. of Eq. ( 13. 3ft . we immediately recognize (after an 
integration by parts) that it vanishes (apart from a total derivative): 



1 



4F1 



d,S 2 x d»S x + p x nS x 



4F X 



0. 



(3.5) 



Therefore, the four-meson Lagrangian (I3.3P reduces to: 



4F2 



+ 



5 

6^2 



Tr [M$ 4 ] 



(3.6) 



In the limit c — Y 0, — )■ and S*x — > this Lagrangian reduces to the usual four-meson 
Lagrangian derived by Di Vecchia et al. in Ref. [23J. 

The last term in the four-meson Lagrangian (I3.6P can be re-written in terms of the mass 
eigenstates, given, in the case A = 0, by Eqs. (I2.27p -( )2.28p . so obtaining: 



5C 



V^ c 1 c 



6 V F W F 



x 



(3.7) 



This term contributes only to the elastic scattering amplitude r] X Tix — > VxVx- At the end 
of the next subsection we shall see that, for A = m u — m d 7^ 0, the term 8C£ gives also 
contributions to the decays into three pseudoscalar mesons, but these contributions are 
strongly suppressed for small A. 
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3.2. The mass eigenstates in the case A/0 



In the strong decays of 77, 77', r]x into three pions the SU(2) isotopic spin / is not conserved, 
i.e. (being the charge conjugation C conserved by strong interactions) the so-called G 
parity, defined, for a multiplet of isotopic spin /, as G = Ce ml2 = C (— l) 7 , C being the 
eigenvalue of C for the neutral component of the multiplet, is not conserved. The mesons 
7], if, r]x are isosinglets (/ = 0) with C = 1 (they can decay into 77 for the electromagnetic 
interaction), and so they have G — 1. On the contrary, the mesons 7r form an isotriplet 
(1 = 1), with Co = 1 (since 7r° can decay into 77 for the electromagnetic interaction), 
and so each of them has G = — 1, and a three-pion final state has G = (— l) 3 = —1. 
We shall evaluate the decay amplitudes (and the corresponding decay widths) at the 
lowest order in the parameter A = m u — m^, which measures isospin violations, i.e., the 
explict breaking of the SU {2) v symmetry. In the case A^O, the fields 7r 3 , ti 8 , S n , Sx mix 
together with the squared mass matrix /C, given by Eq. (I2.24p . while the remaining 7r a 
are already diagonal [8]. We write the matrix K as: 



K = /C + 5/Ca, 



(3.8) 



where /Co is the matrix /C for A = 0, given by Eq. (I2.25p . which is diagonalized by the 
orthogonal matrix Co, given by Eq. (12. 28ft . while 5/Ca is given by: 



/ 



5/Ca 



V 





|SA 




73 5A 







§5A (A 





0/ 



(3.9) 



We shall evaluate the eigenvalues and the eigenstates of the matrix /C at the first order 
in the parameter A, by treating the term 5/Ca as a small perturbation. It is easy to 
verify that the corrections to the eigenvalues (i.e., to the squared masses m^, ml, m 



2 

if • 



m 



evaluated in the previous section) are of order A 2 (the first-order corrections being 
identically zero) and are therefore negligible, if we stop at the first order in A. Instead, 
the eigenstates of the matrix /C at the first order in the parameter A are given by: 







/tt°\ 




(So 


61 


S 2 






= c 


V 


, c = 


a 


Oil 




a 3 




T]' 


A) 


Pi 






\Sx) 




\Vx) 




\7o 


7i 


72 


73/ 



(3.10) 
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where: 



S 
ocq 



A) 



+ 



7o 



+ 



Oil 



a 2 
S 3 
a 3 



1- 

BA 



COS ip 

(ml - m 2 ) 



COS if 



V2F W . 



Frf 



■ Sill ip 



sin (p 



{ml - ml) 



cos ip — sin ip 



BA 



F v sin dp 



{ml - m 2 )Frf 



cos ip + 



F v > 



■ sm <p 



Fn cos <p I y/2F w 



ml - m 2 v ,)F v , y Frf 
Fx sin <p 



cos if — sm ip 



3V2F 2 



x 



(ml — m 2 )F 2 , 



BA 



{ml - m 2 )F v < 



cos ip + 



V2F W . 



■ sm cp 



F x cos I \f2F-x 



(ml - m 2 )F n/ \ F n , 



cos if — sm (p 



V2F w Fx 
(ml-m 2 x )F 2 



BA 



V3(m 2 



- , y/2F n . 
cos tp H — — sm ip 



mi 



Frf 



cos y?, /?! = — sm ip, 71 = — — sm ip, 
frf frf 



BA 



V2F n 



V3(m 2 



mi 



Frf 



cos (p — sm (p 



F n 



-simp, (3 2 = — cos ip, 72 
frf 



V3F 



x 



■ COS ip, 



x 



(m 2 x - ml) Frf 



0, 03 



V3F; 



X 



Frf 



73 



Frf 



(3.11) 



where ml, m\, m 2 ,, m 2 x are given by Eqs. flZZSD , ( I2T29D . (ESQ]), ( IOTP . The only 
modifications with respect to the "unperturbed" matrix C , reported in Eq. (I2.28p . are 
in the elements «o, 0o, 70, £1, S2, S3, which are now different from zero and of order A: 
in the limit A — > the matrix C correctly reduces to the matrix Co- 
At the end of the previous subsection we had observed that in the case A = the last 
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term 8C£ in the four-meson Lagrangian (13. 6p . being proportional to T] x , contributes only 
to the elastic scattering amplitude rjxTjx — > VxVx- Instead, in the realistic case in which 
A = m u — ma 7^ 0, this term has the form (obtained using Eqs. ( 13.10l) -( j3.11l) derived 
above): 



5C 



c / V3 q 1 \ c ( V2BA „ 



-1 4 



VX' 



F^F 



Vx 



x 



(3.12) 



Therefore, when A 7^ this term contributes also to the decay rjx 37r°, but this 
contribution is of order 0(A 3 ), and therefore it is strongly suppressed, for small A, when 
compared with the similar contributions derived from the other terms in the Lagrangian 
(13. 6p [see Eq. ( 13 .211) below]. Therefore, in the following we shall neglect this contribution. 



3.3. Decays 77, 7/, r/x — > 37r°, 



TV TV TV 







In this section we shall evaluate the leading-order (LO) amplitudes and the corresponding 
widths for the decays of rj, rf and t]x into 3tv° or tv + tv~tv°. The fields in the four-meson 
Lagrangian £ 4 , written in Eq. (13. 6p . can be expressed in terms of the fields of the physical 
eigenstates (which diagonalize the squared mass matrix fC) by using Eqs. (13. 10jl — (13.1 ip . 
Let us start considering the decay 77 — > 3tv°. As we have already said after Eq. ( 13. 41) . 
the first term (containing field derivatives) of the four-meson Lagrangian £ 4 in Eq. (13. 6 p 
does not contribute to this decay amplitude, which, therefore, turns out to be simply a 
constant, i.e., not dependent on the particle momenta, and given by, at the first order in 
the parameter A: 



A( V -> 3tt°) = (Tv Tv TV \C A \r]) 



B 



V3F* 



|a(«! + v^i) + 2V3m ^ + («! + V2/3i)(a + ^0)] } • (3-13) 



Using the expressions (I3.1ip for ai, Pi, 5i, ao, 0o and expanding up to the first order in 
the quark masses, we obtain the following expression: 



A(r] -> 37T°) 



BA 
7^ 



cos (p H — — sin ip 



F v ' 



(3.14) 
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The amplitude for the decay r/ — > 3n° can be obtained by simply substituting (Si, ati, (3i) 
with (S 2 , 0(2, P2) m the expression (13.131) . We thus obtain the following expression: 



A(rj' ->• 3tt°) 



BA 

V^F 2 



V2F W 



cos if — sin <p 



(3.15) 



Let us observe that in the limit F x — > (that is, F v > — >■ F w ) the expressions (13. 14j) and 
(I3.15P correctly reduce to the corresponding expressions derived by Di Vecchia et al. in 
Ref. E31 i.e.: 



A(77->37r°)| Fx=0 
A(7/-^37r°)| Fjr=0 



BA 

T^F 2 

BA 



I COS (f + sin if I 
cos if — sin 9? 



(3.16) 
(3.17) 



where (p is the mixing angle without the contribution coming from the new U(l) axial 
condensate, and it is given by Eq. (II. 8p with Fx = 0, i.e., tany? = ^jBF 2 (m 



m) 



F 2 



6V2A^ m v 

From the amplitudes (13.141) and (I3.15p . we can derive the corresponding decay widths by 
integrating over the final-state phase space, according to the formula (valid for constant 
amplitudes A and three identical final particles) F — ^ / ^,d$^ 3 ^|y4| 2 = 2 ^ 3 , $^ 3 \ where 
the total phase space $ ( - 3 ^ is given by (see, for example, Ref. [30J and references therein): 



ml). 



* (3> = / iw = iw r -..)(. -*)(.» -.)(..-.), (3.18) 

where M is the mass of the initial particle and Si = (mi — m 2 ) 2 , S2 = + WI2) 2 , 
s 3 = (M — ms) 2 , S4 = (M + m^) 2 , mi, m 2 and 777,3 being the masses of the three final 
particles; s and t are the usual Mandelstam variables, defined as s = (P — p\) 2 and 
t = (P —P2) 2 , where P is the four-momentum of the initial particle and pi, P2, P3 are the 
four-momenta of the three final particles (P 2 = M 2 , p\ = m\, p\ = m\, p\ = m|). 
After performing numerically the integration in Eq. (13.181) for the two cases that we are 
considering, using the values for the meson masses as reported by the Particle Data Group 
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2 $0) 


$(3) 




m v 


2 $(3) 


$(3) 


rrirf ' 





we have obtained the following expressions for the decay widths: 
r L o(^37r°) = ^cos^ + ^|^sin^ ] ^— , ^— = 9.82 keV, (3.19) 

iW^Svr ) = ^|£^^cos^-sin^] ^- = 67.00 keV(3.20) 

Proceeding analogously, the following expression is obtained for the leading-order ampli- 
tude of the decay rjx — > 37r°: 

-4fc^3,r») = ^^. (3.21) 

7T V 

Let us observe that this amplitude correctly reduces to zero when Fx — > 0, i.e., when 
the new £7(1) axial condensate is zero. Concerning the derivation of the decay width, the 
mass m r)x of the exotic meson rjx is not directly known (but see the discussion at the 
end of this subsection) and therefore the integration in Eq. (13.181) cannot be performed 
numerically. However, on the basis of what we have said in the previous section, the mass 
of the rjx is expected to be quite large, at least larger than the mass of the rf. So it is 
probably not a too bad approximation to neglect the meson masses in the total phase 
space for this process. In the limit m x = m 2 = m 3 = Eq. (13 . 18[) reduces to 

and for the width of the decay rjx —> 37r° we obtain the following approximate expression: 

T LO (Vx 3tt°) = \A( Vx 3n°)\ 2 ^^4 = (^ A f^ m ( 3 . 23 ) 

Let us now study the decays of rj, rj' and rjx into 7r + 7r~7r°. As we have already observed 
above, also the four-meson Lagrangian term 5C^/\ defined in Eq. H3.4j) and contain- 
ing derivatives of the fields, gives contributions to the amplitudes of these decays. In 
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particular, one finds that: 



5A^\r]' ^tv+tv'tv ) 

5A {f \ Vx -> TT + TT-TT ) 



7T 7T 7T 



5£ 



(/) 



TV TV TV 



5C 



(/) 



V 
V 

Vx 



F2 



5 Q 5i(s - s ), 



j^S S 2 {s - s' ), 



) = ^5 5 3 (s-s ), (3.26) 



(3.24) 
(3.25) 



where the coefficients 5 , Si, 5 2 and 5 3 are defined in Eqs. (I3.10j) - (l3.1ip . while s , s' and 
So are so defined: 

s o = ^(m 2 v + 3m 2 n ), s' = i(mj + 3m 2 ), s = -(m 2 ^ + 3m 2 ), (3.27) 

and, as usual, s = (P — P^o) 2 = (P n + + P-^-) 2 , P being the four-momentum of the initial 
particle (rj, rj', rjx) and P n o, P n +, P n - being the four-momenta of the three final pions. 
Adding also the contributions coming from the second term in the r.h.s. of Eq. (I3.6p . we 
obtain the following expressions for the amplitudes of the decays rj, rf, rjx — > tv + tv~tv°: 



A(fj — > TV + TV 7T°) 
A{r)' -> TV + TV~TV°) 
A(rj x -> TV + TV~TV°) 



BA 



COS if 



3a/3F2 

BA I y/21 \ 
3^ V 
V2BAF X 



F v > 



■ sin if 



F n < 



■ cos </? — sin (p 



1 + 



1 + 



3(s - go) 
m 2 — ml 

jKg - ap) 
m 2 ; — m 2 

7?' 7T 



1 + 



3(s - sp) 
ml — ml 

r/x k 



(3.28) 
(3.29) 
(3.30) 



From these amplitudes we can derive the corresponding decay widths by integrating over 
the final-state phase space, according to the formula (see, for example, Ref. [30J and 
references therein): 



2M 



d^\A\ 2 



dsdt 



2M / 128vr 3 M 2 



\A\ 



I f S3 ds 

256tt~W J, -|^)lV(s-si)(s-s 2 )(s 3 -s)(s 4 -s), 



(3.31) 



where the notation is the same already used in Eq. (j3.18p . After performing numerically 
the integration in Eq. (I3.3ip . using the values for the meson masses as reported by the 
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Particle Data Group [21], we have obtained the following expressions for the decay widths: 
Tlo(v -> vt+tt-tt ) = Lo S ^ + ^!L s in^ x 10.48 keV, (3.32) 

r L o(V -> 7r + 7r"7r ) = ^f^l f cosy -sin y I x 83.95 keV. (3.33) 

Concerning the case of the decay rjx — > 7r + 7r~7r°, we proceed exactly as for the case of 
the decay rjx — > 37r° and we neglect the meson masses in the calculation of the integral 
(I3.3ip . so obtaining the following approximate expression for the decay width: 



T LO (vx vr+TT-vr ) = , i ^ m,,. (3.34) 



15367r 3 FiF 2 ,"^' 



7r r\' 



We now numerically compute our theoretical expressions for the leading-order decay 
widths, using for the mixing angle (p the value derived from Eq. (jl.8p . 
All our isospin-violating decay widths are proportional to the factor: 



(BA) 2 = m A ( m " m<i ) 2 = mi ( ^-^) " ~ 2.66 x 10 7 MeV 4 , (3.35) 
7T \m u + m d J n \R+lJ 

where m\ = B(m u + ma) — (134.98 MeV) 2 and R = m n /m^ ~ 0.558 is the ratio between 
the up and down quark masses, determined using Eqs. (I2.23P and the experimental values 
of the meson masses reported in the Particle Data Group [24]. 

We are, of course, particularly interested in the effects due to a non-zero value of the 
parameter Fx, related to the new U(l) axial condensate considered in this paper. In Table 
1 we report, for each decay process of the form 77,77' — > 3tt, the leading-order theoretical 
prediction, using for the parameter F x the value F x = 24(7) MeV, that we have found 
studying the radiative decays r], 7/ — > 77 [see the Introduction and, in particular, Eq. 
(11.91) ]. These values are compared with the corresponding values obtained for F x = 0, 
i.e., in the absence of the new U(l) axial condensate (in Table 1 we also explicitly show 
the correction to the decay widths, AFlo = ^ho(Fx = 24 ± 7 MeV) — IYo^x = 0), 
coming from a non-zero value of Fx), and also with the experimental values. 

Concerning the comparison with the experimental values, it is well known that, because 
of large unitarity corrections due to strong final-state interactions, one has to go beyond 
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Dpchv 


exp V v / 


F x = F x = 24(7) MeV 


Arm (keV) 

i—ij. ^ivc; v y 


7] -> 3tt° 


0.423(26) 


0.178 0.176(1) 


-0.002(1) 


7/ ->■ 3tt° 


0.33(6) 


0.84 0.62(10) 


-0.24(10) 


7] — )■ 7T + 7r~7r° 


0.30(2) 


0.127 0.125(1) 


-0.002(1) 


77' — )■ 7T + 7r _ 7r° 


0.70(25) 


0.70 0.52(8) 


-0.18(8) 



Table 1: The leading-order theoretical predictions for the decay widths, computed both 
for Fx = and for Fx = 24(7) MeV, and the corresponding corrections to the decay 
widths, ArYo = Tlo(Fx — 24 ±7 MeV) — Tj J o(Fx = 0), compared with the experimental 
values. 

leading and even one-loop order in chiral perturbation theory in order to obtain a valid, 
reliable representation of the 77, 77' — > 3n decay amplitudes and of the corresponding decay 
widths, that can be successfully compared with the experimental values [3H [32j [331 EH 
[351136]. 

In the present paper we are not, of course, aiming at that. In particular, we cannot 
proceed as in the case of the radiative decays 77,77' —> 77, i.e., we cannot extract the 
value of Fx (and of the mixing angle (p) by comparing, e.g., the leading-order theoretical 
predictions (13. 19ft and (I3.20p . for the rj — > 3tt° and 7/ — > 3tt° decay widths, with the 
corresponding experimental values reported in Table 1. Indeed, making use of Eq. ( II. 8p 
for tan<^5, one easily verifies that, being tan <p and F v > = \JF% + 3F| increasine; functions 
of Fx, the expression ( 13.191) for Flo(^ - > 37r°) is a decreasing functions of Fx- so, being 
its value at Fx = already smaller than the corresponding experimental value, it turns 
out that there is no value of Fx which makes the expression ( 13 . 191) compatible with the 
experimental value in Table ljj 

Instead, our aim is simply to quantify the corrections coming from a non-zero value of 
the parameter F x , taking the leading-order amplitudes/widths in the F x = case as a 
useful reference point. From the values reported in Table 1 we can conclude that: 

i) In the case of the 77 — > 3tc decays, the size of the corrections AI^lo coming from a 

*Even considering the singlet decay constant F v > and the mixing angle (p in Eqs. p. 191) and (|3.20[) 
as free parameters, to be fixed from a comparison with the experimental values reported in Table 1, we 
would find a too small value F n * ~ 68 MeV for the singlet decay constant, incompatible with the formula 
(TO]) , i.e., F v i = ^fFj+~3Fj > Fyr = 92.2(4) MeV, and also an anomalously large value (p ~ 44° for the 
mixing angle. 
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non-zero value Fx = 24(7) MeV, with respect to the Fx = case, is very small, 
being of the order of 1%, i.e., comparable to (or even smaller than) the size of the 
electromagnetic corrections for these decays, which have been recently re-calculated 
inRef. [37]. 

ii) Instead, in the case of the rf — > 3n decays, the size of the corrections Ar L o is 
much larger, being of the order of 30%. Moreover, at least for the decay rj' — > 37r° 
(where the statistical errors are smaller), this (negative) correction seems to go in 
the right direction, improving the agreement between the theoretical prediction and 
the experimental value. 

Concerning the decays of the rjx into three pions, we derive the following relations 
between its mass m rix and the decay widths: 

rLo( ^ 37r ° } = Flo( ^ ~* n+ " 0) = (4.35i?#) x 10- 7 . (3.36) 

These constraints could be used to identify a possible candidate for the exotic singlet 
meson rjx, once we know its mass and decay widths. According to the Particle Data Group 
|24j . the possible candidates for the rjx, having the same quantum numbers I G (J PC ) = 
+ (CT + ) of the rj', but a larger mass, are the following: 



77(1295) 


"tot 


= 55(5) MeV, 


77(1405) 


Ttot 


= 51(3) MeV, 


77(1475) 


Ttot 


= 85(9) MeV, 


77(1760) 


Ttot 


= 96(70) MeV, 


77(2225) 


Itot 


= 18512 MeV - 



(3.37) 

Unfortunately, no quantitative determination of their decay widths into three pions has 
been done up to now. 



3.4. Decays 77' — > i]tttt and rjx — > tjtitt, tj'tttt 

We now study the decays of 77' into 777r 7r°, i]7i + 7i~ and of 77^ into r?7T 7r°, r?7r + 7r~ , 7]'ir 7r , 
ri'7T + 7T~ . This decays do not violate isospin and so they can happen also when A = 0. 
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Therefore, in order to evaluate the amplitudes and the corresponding widths for these 
decays, we shall use the approximate expressions ( I2.28P of the eigenstates at the order 
zero in the isospin-violating parameter A. 

The following expression is obtained for the leading-order amplitudes of the decays rf — > 
7]7i°7T and rj — > i]7i + 7i^ (which, in the limit of exact SU(2)y isospin symmetry, are equal): 



A(j]' r]7r"7r 
,2 



A{rf — > 7/7T 7T~ 



mt 



QF2 



2x/I/; cos(2^)+ f^J- I )sn,(2^) 



F 2 



(3.38) 



In the limit Fx — > these amplitudes reduce to the expression already found in Ref. [23] . 



i.e.: 



2 r 



6F2 



2^008(2^) + sin(2^)j . (3.39) 



After numerical integration of the phace space f 1 3 . X 8 [) . using the values for the meson 
masses reported in Ref. [21], we obtain the corresponding decay widths: 



0\\ 2 



\A{r)' rjTrV) 



$(3) 



$(3) 



2rry • 2! ' 2rry ■ 2! 

9 $(3) 

L4(V -> r/TT+vr-) = 2r LO (r/ -> tjTrV). 

1 1 Im^i 



1.093 keV, 



(3.40) 



We proceed as in the previous subsection and numerically compute our theoretical expres- 
sions for the leading-order decay widths, using for the mixing angle (p the value derived 
from Eq. (jl.8p and for the parameter Fx the value Fx = 24(7) MeV, that we have 
found studying the radiative decays rj, rf — > 77. Again, our aim is simply to quantify the 
corrections coming from a non-zero value of the parameter Fx, taking the leading-order 
amplitudes/widths in the Fx = case as a reference point. In this case, however, it is 
already known from Ref. [23] that the leading-order theoretical predictions for Fx = 0, 



r L o(?/ ->• Vk + k-)\f x =o = 2T w (rf -> r/7r°7r°)| Fx=0 = 2.42 keV 



(3.41) 



are in strong disagreement with the experimental values (24], T exp (r]' — > i]7i + 7i ) = 84(5) 
keV and r cxp (r/ -> r]n°n ) = 42(4) keV. 

We can try to see if the introduction of a non-zero value of Fx can cure, at least in part, 
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the strong disagreement between leading-order theoretical predictions and experimental 
values: however, the answer to this question is negative. In fact, we find that: 



Tlo(V -> r ? 7r+7r )|f x =24(7) Mev = 2T LO (rj' -»■ r^vr ) |f x =24(7) MeV = 1.78(30) keV. (3.42) 

Even if the correction AlTix) is quite large (of the order of 30%) if compared with the value 
of Tlq at Fx = 0, it is, however, too small if compared with the experimental value. In 
addition, the correction AI^lo, being negative, goes in the "wrong" direction, lowering the 
theoretical prediction at Fx = 0, which is already much smaller than the experimental 
value: in other words, it is not possible to find a value of the parameter Fx which moves 
the leading-order theoretical prediction towards the experimental value. Moreover, the 
amplitude (I3.38P is a constant, while the experimental data are well fitted by a non- 
constant amplitude having the form: A(j]' — > rjirir) = A(l — <JiT v ), where T v is the kinetic 
energy of the r), A and <j\ are some constants. As already observed in Ref. [23], in 
order to describe this behaviour, and to obtain a better agreement with the experimental 
value of the decay width, it is not enough to retain only the leading order in the 1/N 
expansion, but one has to go to the next-to-leading order, adding to the Lagrangian ( 11.211 
non-leading terms such as XQ 2 Trld^Ud^W), that may be very important because of the 
proportionality of the leading terms to the tiny pion masslll The systematic introduction, 
in our model, of higher-order terms in the 1/N expansion (including also one-loop graphs, 
which are of order 1/N 2 : see, for example, Refs. [311 [39] ) is, of course, a quite hard task, 
which is beyond the aim of the present paper (but it will probably be addressed in a 
subsequent work). 

Concerning the exotic meson rjx, the following expressions are obtained for the leading- 
order amplitudes of the decays rjx — > V 7171 an d Vx — > rfnn: 

A(r) X 777TV) = A{r lx ^r 1 ^n-)= ^f F X i Lsy+^sinyj , (3.43) 

Ainx^rjW) = A( Vx -> 7W-) = ^f 2 F p X cos^ - siny^ . (3.44) 

From these amplitudes we can obtain the corresponding decay widths, using, for the 
integrated phase space (I3.18p . the following approximate expression obtained neglecting 



different and alternative approach, first suggested in Ref. [35], considers the decay rj — > rjirir to be 
dominated by coupling to nearby scalar resonances. 
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the pion masses (while retaining the rj and rf masses different from zero): 



$j 3) (M,m) 



M 4 - m 4 + 4M 2 m 2 ln(m/M) 
256tt 3 M 2 ' 



(3.45) 



where M is the mass of the initial particle and m is the mass of the final massive particle. 
We thus find the following expressions: 
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(3.46) 



(3.47) 



As in the case of Eq. f)3.36p . also these relations could in principle be used to identify 
a possible candidate for the exotic singlet meson rjx, once we know its mass and decay 
widths. However, a certain caution must be used since, as in the case of the decays 
7]' — > r]7T7T, large corrections to these leading-order results could come from non-leading 
terms in the 1/N expansion: only a detailed analysis of our model at the next-to-leading 
order in 1/N shall clarify this point. 
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3.5. Possible decays r)x 37/, 7/7/7/, rjrj'rf, 3r/ ? 

If the exotic singlet meson r)x were heavy enough, let us say, if m Vx > 3m v ~ 1640 
MeV, it could also decay into three 77 particles. The amplitude for this decay, which does 
not violate SU(2) V isospin, can be evaluated at the order zero in the isospin-violating 
parameter A, so using the approximate form ( I2.28|) for the physical eigenstates, and the 
following result is obtained: 

Mm 3 "> = {~ costp+ -^7 am ' p j ' (3 ' 48) 

In order to estimate the decay width, considering that also the final-state particles 77 are 
rather heavy, we use the approximate expression for the total phase space ( 13.181) in the 
non-relativistic limit, i.e.: 

^(M,m 1 ,m 2 ,m 3 ) = -g^ J (3 _ 4g) 

647T Z y (mi + m 2 + ma) 13 

where Q = M — m\ — mi — m 3 is the so-called U Q value" of the decay. 
We thus obtain the following approximate expression for the decay width: 

C(m, x) m„m„m,) 



T L o(vx = \A(t} X ^3t/) I 2 



2m w ■ 3! 



m%F x ( _ 3v/2^ . (m w -3m,) 2 

1 cos if — sm if 1 



m vx 



486x/3vr 2 F4F2 ^ r 2F V 

(0-96±8;«) x 10" V* f 1 - 2^ V , (3.r )0 ) 



m Vx 



where, as usual, we have used for the parameter Fx the value Fx = 24(7) MeV, that 
we have found studying the radiative decays 77, 7/ — > 77, and for the mixing angle <p the 
value derived from Eq. (|1.8|) . (For example, for a value m Vx ~ 2 GeV, one would get 
r LO fe -> 3t/) « 61 keV.) 

Other possible decays of this kind (supposing that the ryx is heavy enough so that they are 
kinematically allowed) are rjx VVV'i w'v'i 37/', and their amplitudes and corresponding 
widths can be derived in a similar way. 
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4. Conclusions 



In this paper we have considered a scenario (supported by some lattice results) in which 
a £/(l)-breaking condensate survives across the chiral transition at T c h, staying differ- 
ent from zero up to Tjj^ > T ch . This scenario has important consequences for the 
pseudoscalar-meson sector, which can be studied using an effective Lagrangian model, 
including also the new £7(1) chiral condensate. This model, originally proposed in Refs. 
[SI [71 [H [9] and elaborated in Refs. [TUl [III [12], could perhaps be verified in the near future 
by heavy-ion experiments, by analysing the pseudoscalar-meson spectrum in the singlet 
sector. 

Section 2 contains a brief review (for the benefit of the reader) of the main results, 
obtained in the original papers [61 EJ [10], concerning the mass spectrum of the Chiral 
Effective Lagrangian. The Lagrangian (12 .ip contains a new field X and three new pa- 
rameters, namely Fx, 0J\ and ci, with respect to the usual Lagrangian of Witten, Di 
Vecchia, Veneziano et al. In this paper we have assumed that the parameter Fx, which 
is essentially proportional to the new U(l) axial condensate, is different from zero. In 
this case, there are two singlet pseudoscalar mesons, the rf and an exotic particle rjx, 
whose squared masses (assuming also that the coupling constant c\ of the interaction 
term det(L r )X"'" + det(U^)X in Eq. ( 11.31) . between the usual qq meson field U and the 
exotic meson field X, is different from zero and not too small: see the discussion in Ap- 
pendix B) are given by Eqs. (I2.30p and (I2.3ip : in particular, the exotic particle rj X turns 
out to have a large (non-chiral) mass term of order 0(1) in the large- N limit, generated 
by the (non-zero) coupling constant C\. 

In section 3, generalizing the results obtained in Refs. [TH [12], where the effects of 
the new £7(1) chiral condensate on the radiative decays of the pseudoscalar mesons r\ and 
rj' into two photons had been investigated, we have studied the effects of the U{1) chiral 
condensate on the strong decays of the "light" pseudoscalar mesons, i.e., rj,rj' — > 3n°; 
T],rf — > 7r + 7r~7T°; r)' — > i]7c°tt ] rf — > i]'K + 7i^; and also on the strong decays of the exotic 
("heavy") 5*?7(3)-singlet pseudoscalar state rjx- r\x —> 37r°; r/x — > 7r + 7r~7r ; rjx —> ^7r°7r°; 
f]x rj7i + 'K~; r]x —> rj'^n ; r/x —> 7/7r + 7r~; rjx —> 3r/, 3r/, r/r/r]', r/r/'r]'. Concerning the 
decays of the exotic particle r]x, we have found some relations between its mass and its 
decay widths, which in principle might be useful to identify a possible candidate for this 
particle. According to the Particle Data Group [24] . the possible candidates for the rjx, 
having the same quantum numbers I G (J PC ) = + (0~ + ) of the rj', but a larger mass, are, 
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at the moment, those reported in Eq. (13.371) (other candidates with larger masses are also 
present, but some of their quantum numbers I G (J PC ) are not yet known): unfortunately, 
no quantitative determination of their decay widths into (e.g.) three pions has been done 
up to now. 

Concerning the decays T],T]' — > 3ti, it is well known that, because of large unitarity 
corrections due to strong final-state interactions, one has to go beyond leading and even 
one-loop order in chiral perturbation theory in order to obtain a valid, reliable repre- 
sentation of the decay amplitudes and of the corresponding decay widths, that can be 
successfully compared with the experimental values [3TJ [32j [33J, [3U [351 ES] . 
In the present paper we have not, of course, aimed at that. In particular, we could not 
proceed as in the case of the radiative decays 77,7/ — > 77, i.e., we could not extract the 
value of Fx (and of the mixing angle (p) by comparing, e.g., the leading-order theoretical 
predictions ( 13 . 1 9[) and (I3.20p . for the 77 — > 3n° and 7/ — > 37r° decay widths, with the 
corresponding experimental values reported in Table 1 of section 3. 
Instead, our aim has been simply to quantify the corrections coming from the non-zero 
value Fx = 24(7) MeV, that we have found studying the radiative decays 77, 7/ — > 77 [see 
the Introduction and in particular Eq. (II .9)1 ]. taking the leading-order amplitudes/widths 
in the Fx = useful reference point. From the values reported in Table 1 of 

section 3 we have concluded that: 

i) In the case of the rj — > 3tc decays, the size of the corrections AI^lo coming from a 
non-zero value Fx = 24(7) MeV, with respect to the Fx = case, is very small, 
being of the order of 1%, i.e., comparable to (or even smaller than) the size of the 
electromagnetic corrections for these decays, which have been recently re-calculated 
in Ref. [37]. 

ii) Instead, in the case of the 77' — > 3tt decays, the size of the corrections AI^lo is 
much larger, being of the order of 30%. Moreover, at least for the decay r/' — > 3n° 
(where the statistical errors are smaller), this (negative) correction seems to go in 
the right direction, improving the agreement between the theoretical prediction and 
the experimental value. 

Finally, concerning the decays 7/ — > 777r°7r° and 7/ — > r/7r + 7r _ , knowing already from Ref. 
[23] that the leading-order theoretical predictions for Fx = are in strong disagreement 
with the experimental values, we have tried to see if the introduction of a non-zero value 
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of Fx can cure, at least in part, this disagreement: but we have found that it cannot. In 
fact, even if the correction Ar L o is quite large (of the order of 30%) if compared with the 
value of r L o at F x = 0, it is, however, too small if compared with the experimental value, 
and, moreover, being negative, it goes in the "wrong" direction, lowering the theoretical 
prediction at Fx = 0, which is already much smaller than the experimental value. (In 
other words, it is not possible to find a value of the parameter Fx which moves the 
leading-order theoretical prediction towards the experimental value.) 

However, as we have already stressed in the conclusions of Refs. [HI [T2], one should 
keep in mind that our results have been derived from a very simplified model, obtained 
by doing a first-order expansion in 1/N and in the quark masses. We expect that such 
a model can furnish only qualitative or, at most, "semi-quantitative" predictions. As 
already observed in Ref . [23] , in order to obtain a better agreement with the experimental 
data of the decay widths, most probably it is not enough to retain only the leading order 
in the 1/N expansion, but one has to go to the next-to-leading order. The introduction, 
in our model, of higher-order terms in the 1/N expansion is, of course, a quite hard task, 
which is beyond the aim of the present paper. Further studies are therefore necessary in 
order to continue this analysis from a more quantitative point of view. We expect that 
some progress will be made along this line in the near future. 
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Appendix A: The U(l) chiral order parameter 



We make the assumption (discussed in the Introduction) that the U(l) chiral symmetry 
is broken independently from the SU(L) ® SU(L) symmetry. The usual chiral order 
parameter (qq) is an order parameter both for SU(L) £g> SU(L) and for U{1)a- when 
it is different from zero, SU(L) <g> SU(L) is broken down to SU(L)y ("V stands for 
"vectorial") and also U(1)a is broken. In fact, under a U(l) chiral transformation with 
parameter a (as usual, qi = |(1 + 75)9 and q R = |(1 — 75)g, with 75 = — Z7°7 1 7 2 7 3 , 
denote respectively the left-handed and the right-handed quark fields): 

U(1) A : q^e- ia ^q, i.e., q L ^ e~ ia q L , q R -> e ia q R , (A.l) 

the chiral condensate would transform as (assuming the U(1)a symmetry to be realized 
a la Wigner-Weyl) : 

U(1) A : (qq) e 2 ^(q L q R ) + e~ 2ia (q R q L ). (A.2) 

By taking a = tt/2, we would obtain (qq) — > —(qq): therefore, if the chiral condensate 
is different from zero, the U(1)a symmetry cannot be realized a la Wigner-Weyl. Thus 
we need another quantity which could be an order parameter only for the U(l) chiral 
symmetry [6j [71 EJ |9] . The most simple quantity of this kind was introduced by Kobayashi 
and Maskawa in 1970 [16] . as an additional effective vertex in a generalized Nambu-Jona- 
Lasinio model, and it was later derived by 'tHooft in 1976 [3], as an instanton-induced 
quark interaction. (See also Ref. [17] for an historical review on this subject.) 
For a theory with L light quark flavours (of mass rrii Aqcd] i = 1, • • • ,L), it is a 
2L-fermion interaction that has the chiral transformation properties of: 

+ H.c. = det (q sR q tL ) + det (q sL qtR) , (A.3) 

st st 

where s,t = 1, . . . ,L are flavour indices, but the colour indices are arranged in a more 
general way (see below). Since under chiral U(L) ®U(L) transformations the quark fields 
transform as follows: 

U(L)®U(L): q L ^V L q L , q R ^ V R q R , (A.4) 
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U U(1) 



det 

st 
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where Vj, and V R are arbitrary L x L unitary matrices, we immediately derive the trans- 
formation property of under U(L) <8> U(L): 

U(L) ® C/(L) : ©g^ -> det(Vb) det(Vfi)*det (q sR q tL ) + H.c. (A.5) 

This just means that OEL is invariant under SU(L) <8> SU(L) <g) [7(l)y, while it is not 
invariant under the {7(1)a transformation (lA.lj) : 

C/(1) A : -> e" i2La det (g si? g iL ) + H.c. (A.6) 

A.l The U(l) chiral condensate for L = 2 

As an example let us consider the most simple case, that is L = 2, but with a general 
colour group SU(N). It is not hard to find (using the Fierz relations both for the spinorial 
matrices and the SU(N) generators in their fundamental representation) that the most 
general colour-singlet, Hermitian and P-invariant local quantity (without derivatives) 
which has the required chiral transformation properties is just the following four-fermion 
local operator: 

O^f K A)) = K 0o)e st {q a 1R q b sL • q c 2R qf L + q a 1L q b sR • q c 2L qf R ) , (A.7) 

where the colour tensor F b a d (a , (3 ) is given by: 

FSftao, ft) = a S a b 6 d + /3 «, (A.8) 

«o and /?o being arbitrary real parameters. In Eq. (1A.7[) . a,b,c,d G {l,...,iV} are 
colour indices; s,t e {1,2} are flavour indices and e s * = — e ts , e 12 = 1. Dirac indices are 
contracted between the first and the second fermion field and also between the third and 
the fourth one. Note that if we choose ao = N and ft = —1, OL^ («o, ft) just becomes 
(up to a proportionality constant) the effective Lagrangian for two flavours of quarks in 
an instanton background, found by 'tHooft in [3] . 

Now, to obtain an order parameter for the U(l) chiral symmetry, one can simply take 
the vacuum expectation value of O v ^ (a , ft): 

fi'K, ft) = (OgSfW ft)}. (A.9) 
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The arbitrarity in the choice of «o an d Po (indeed of only one of them, since only their 
ratio is relevant) can be removed if we require that the new U(l) chiral condensate is 
"independent", in a sense which will be explained below, of the usual chiral condensate 
(qq). As it was pointed out by Shifman, Vainshtein and Zakharov in [ID], a matrix element 
of the form (qTiq ■ gl^g) has, in general, a contribution proportional to the square of the 
vacuum expectation value (v.e.v.) of qq. This contribution corresponds to retaining the 
vacuum intermediate state in all the channels and neglecting the contributions of all 
the other states; we call this contribution the "disconnected part" of the original matrix 
element: 

(qT iq ■ qT 2 q) disc = -L [(Trl^ • Tr T 2 ) - Tr^)] (qq) 2 , (A.10) 
where the normalization factor G is defined as (qq = J2a^ a( 1 a )- 

(UIb) = -^-{qq), i-e., G = 5 A a, (A.H) 

and the subscripts A, B are collective indices which include spin, colour and flavour; 
therefore, G = 4 x L x N for a general L, and G = 8N for L = 2. When considering the 
operator O^^(a , /3 ) defined in Eqs. f]A.7|) and (1A.8[) . we find the following expression 
for its disconnected part: 

(O[f= 2) (ao, 0o)) di sc = ^l N ( 2a o + A)) + («o + 2/3 )](gg) 2 , (A.12) 

where: (qq) = (uu) + (dd). From this last equation we immediately see that the discon- 
nected part of the condensate («o> Po) vanishes with the following particular choice 

of the coefficients ao and /?o (only their ratio is really relevant): 

a N + 2 y ' 



In other words, the condensate (1A.9|) with a and (3 satisfying the constraint (1A.13j) does 



not take contributions from the usual chiral condensate (qq). To summarize, a good choice 
for a U(l) chiral condensate which is really "independent" of the usual chiral condensate 
(qq) is the following one (apart from an irrelevant multiplicative constant): 

C { u { f = ((m - ^±^«e st (fmQsL • &r4l + qUsR • (A.14) 
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As a remark, we observe that the condensate C^?* so defined turns out to be of order 
0(g 2 N 2 ) = O(N) in the large- N expansion (this result was also derived in Ref. [7] by 
simply requiring that the 1/N expansion of the relevant QCD Ward Identities remains 
well defined when including this new condensate). In the case of physical interest, i.e., 
N = 3, the condensate Eq. ( 1A.14I) becomes: 

4« 2) = am - \my l • + r 1L q b sR ■ (A.15) 

A. 2 The U(l) chiral condensate for L = 3 

So far we have considered the most simple case L = 2. However, this procedure can be 
easily generalized to every L, and we can take as an order parameter for the £7(1) chiral 
symmetry: 

<#(!) = (<!))• (A.16) 
As we have done in the case L = 2, the colour indices may be arranged in such a way that 
the U(l) chiral condensate does not take contributions from the usual chiral condensate 
(qq): as a consenquence of this, the new condensate will be of order 0(g 2L ~ 2 N L ) = O(N) 
in the large- N expansion [8] . 

In the real-world case there are L = 3 light flavours, u, d, and s, with masses m u , rrid, 
and m s which are small compared to the QCD mass-scale Aqcd- Proceeding as in the 
case L = 2 [see Eq. (I A. 71) ]. one reduces to consider the following general colour-singlet, 
Hermitian and P-invariant local six-fermion operator (without derivatives): 

WL=3) _ paia 2 a 3 hhh ^ 01 „*>i . ^ a 2 b 2 - a 3 b 3 , n ( A 17) 

u u(i) ~ r b ± b 2 b z e qm qi X L ^2/? Qi 2 l qm y«3i + rLX -' k^- 1 ') 

where a±, 02, a 3 , b\, 62, 03 e {1, 2, • • • N} are colour indices, h, 12, h € {1, 2, 3} are flavour 
indices and the colour tensor F^^ 3 is given by: 

+ + + w» (a. is) 

with ai, a2, 0:3, /3i, /?2, /?3 real parameters. However, differently from the case L = 2, the 
operator O^ 1 ^ in Eqs. (1A.17j) - (1A.18j) . with arbitrary real parameters ai, a 2 , a 3 , Pi, f3 2 , (3 3 , 
is not, in general, invariant under a (3) <g> SU(3) chiral transformation: 

SU(3)®SU(3): q L ^U L q L , q R ^U R q R (A.19) 
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(det Ul = detUn = 1). Invariance under SU(3) (g) SU(3) is, instead, recovered provided 
that the colour tensor F^ 2 \® 3 satisfies the following symmetry property: 

KZT = F &XC> V Permutations k} of {1, 2, 3}. (A.20) 



In fact, in this case it is easy to see that the operator flA.17j) can be re-written in the 
following form: 

<T) 3) = KZT ±e r ™e lM ° g r - ■ qj» q% L ■ q r f R ft L + H.c, (A.21) 
which is manifestly invariant under £77(3) ® £77(3): 

^n(L=S) paia 2 a 3 ^ rir 2 r- 3 Jihh?; ai (Tj\ \ (TTA, n bl 

W U(X) ~^ r b 1 b 2 b z g| e e ( is 1 R\ U R)s 1 r 1 \ u L)l 1 m 1 q miL 

■ Qs a 2 UU R )s 2 r 2 (U L ) hm2 qt 2L ■ q s ?R(rf)s 3 rs(U L )l 3 m 3 <& L + H.C. 

= K™^ det K) eSlS2S3 det(f/ L )e mim2m3 



x 



i ls 1 R i lrn 1 L "s 2 R Hm 2 L i ls 3 R < dm 3 L ' n ' L - 

T^aia 2 a 3 1 s 1 s 2 s 3 m 1 m 2 m 3 - a ± b 1 - a 2 b 2 - a 3 b 3 , tt 
r b 1 b 2 b 3 01 e e ( is 1 R ( dm 1 L ' 1s 2 R VrnaL ' y S3 -R Vm 3 L n c - 



3! 

) 



^m 3) - (A.22) 



(Or, equivalently, one can start from the expression (IA.21I) of the six-fermion local oper- 
ator, which, on the basis of (IA.22I) . is invariant under SU(3) <S> SU(3) for every choice of 
the colour tensor F^^\ but one immediately recognizes that only the symmetric part 
of the colour tensor, satisfying the relation (IA.20[) . contributes to the r.h.s. of flA.21[) . the 
anti- symmetric parts being trivially cancelled out. Note that, in the case L = 2, the most 
general colour tensor (1A.8I) automatically satifies the symmetry property, = F^.) 
The symmetry property ()A.20j) imposes the following constraints on the parameters of 
the colour tensor flA.18j) : 

« 3 = « 2 , p 3 = fo = p 1 . (A.23) 
The colour tensor has, therefore, the following form: 

+ mtlWl+Wt^t + KlWl), (A.24) 



36 



in terms of three arbitrary real parameters cui, a 2 , 0\ 



Let us now evaluate the vacuum expectation value of the operator O^,^ 



C^f = (<7) 3) ) = KZ> hhh (ffTi? • ?r 2 g ■ gr 3 g) + cc. 

fel&2&3 



F 6 Xt 3 e' 1 ' 2 ' 3 (r 1 )^(r 2 ) Ci3 (r3) Ej F(g J 4gBgc?Dfe^) + cc, (A.25) 



where [see Eq. (1A.17I) ]: 
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® (<W<Wi) ® {6 ciai 5 dlbl ), 

® (<W<W 3 ) ® (A.26) 



where 2,j are Dirac indices, m,n are flavour indices, and c,d are colour indices. 

As in the case L = 2 treated above, we can write the vacuum expectation value of the 

(L=3) 

operator Ojjn\ as the sum of a connected part, which does not depend on the chiral 
condensate (qq), and a disconnected part, which instead contains the chiral condensate 
(qq), i.e., C^ 3) = (0^ ] ) conn + (0 { ^) disc , where: 

(O^)^ = F fe XC eili2/3 ( r i)^( r 2)^(r 3 )EF(gAgBgctog£;gF)d iSC + cc, (A.27) 

and the disconnected part of the v.e.v. of the six-fermion operator has the following form: 

{qAqBqcqDqEqF)disc 

= (qAqB)(qcqDqEqF) conn + (qcqD) (qAqsqEqF) conn + (qEqF)(qAqBqcqn) conn 

-(qAqD)(qcqBqEqF)conn ~ (<M?f) (^D^b) conn ~ (<?C9s) (^D^f) conn 
-(qcqF)(qAqBqEqD)conn ~ (^b) (qAqFqcqD) conn ~ (q~EqD) (^B^f) conn 

+ (qAqB){qcqD}{qEqF) - (uQb) (qcqF) (qEqD) ~ (9a9d> (<?cgs) (^f) 

+{qAqo){qcqF){qEqB) + (uqF) (<7c<?b) (feto) - (uqF)(qcqD)(qEqB)- (A.28) 



On the basis of Eq. ( 1A.11I) . we see that the disconnected part of the condensate (1A.25I) 
can be written as: 

(0$3°>*c = A!<gg> + A 3 (qq) 3 , (A.29) 
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where the first term (Ai(qq)), proportional to the chiral condensate, is originated by the 
first nine terms in the r.h.s. of Eq. (IA.28j) . while the second term (A 3 (qq) 3 ) is originated 
by the last six terms in the r.h.s. of Eq. flA.28j) and represents the completely disconnected 



part, proportional to the third power of the chiral condensate. 

Explicitly, using Eq. i rA~TT]) . with G = 4x3xN = 12N, the form flA~26|) of the T matrices 



and the form flA.24p of the colour tensor Fj££^ s , satisfying the symmetry property (IA.20|) . 
we obtain the following expression for the coefficient A\\ 

conn 

+ Tr T 3 (qT 1 q ■ qT 2 q) 

conn 

- (qT 1 T 2 q ■ qT 3 q) 

conn 

(qT 2 T 1 q ■ qT 3 q) 

conn 

- (qT^rt ■ qT 2 q) conn - (qT^^ ■ qT 2 q) 

conn 

- (qT 2 T 3 q ■ qT^conn - (qT 3 T 2 q ■ qT^)^} + c.c. 
i 

To i^ace i ryace , rjiacei I /~iabcd , /~iabcd , /^iabcd\ ( a on\ 

V' t bde + fbed + fedbl l C 12 + C 13 + G 23 J> (A.3UJ 



12N 



where: 



C" 2 = e st3 [(qiR q h sL ■ q 2R qf^conn + c.c] , 

(jabcd = e s2t j-^-^o ^ . ^ q tL )conn + c - c -] > 

C abcd = e lst [(q 2R q b sL .q 3R qf L ) conn + c.c.], (A.31) 
and the following expression for the coefficient A 3 \ 

As = J^K^r^'^tTrriTrrsTrra-TrriTr^^-Tr^r^Trrs 
+ r Lr(r 1 r 3 r 2 ) + Tr^r^) - Tr^) Tr r 2 ] 

1 - [«! {2N 3 + 3N 2 + N) + a 2 (N 3 + 6N 2 + 5N) 



216AT 3 1 

+ ^(3iV 3 + 9iV 2 + 6iV)]. (A.32) 

Now, if we want to obtain a new order parameter which is really independent on the 
usual chiral condensate (qq), we must require that its disconnected part (IA.29I) vanishes 
independently on the value of (qq), imposing the two conditions A\ = and A 3 = 0. 
Therefore, we have two independent constraints on the three parameters ai, ct 2 and 
0i, which enter the colour tensor ( 1 A. 241) : the new condensate is then univocally 

determined, apart from a multiplicative constant. 
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Let us also observe that in the large- N limit, taking the coefficients a±, a2 and f3\ in the 
colour tensor (IA.24j) to be of order O(N ), the coefficient A\ is of order O(N), while the 
coefficient A 3 is of order O(N ): and, consequently, the first term Ai(qq) in the r.h.s. 
of (1A.29I) is of order 0(N 2 ), while the second term A^(qq) 3 is of order 0(N 3 ) (being 
(qq) = O(N)). If both these disconnected parts are zero, then the new condensate 
is simply equal to the connected part (0\j^ ) C onn, which is of order 0(N), i.e., of the 
same order of the usual chiral condensate (qq) (as already observed in Refs. [7J [8]). 
We also observe that the condition Ax = implies that the new six-fermion condensate 
does not take contributions from four-fermion condensates of the form (IA.31I) . 
In this paper we have only studied the effects of the new six-fermion U(l) chiral order 
parameter. However, recently, four-fermion operators (which could be associated with the 
above-mentioned four-fermion condensates) have been used in the literature, in the study 
of scalar mesons, which are modelled as four-quark (i.e., qqqq) states, called "tetraquarks" 
or " diquark-antidiquark" bound states [1T| W2 \ H3"]. 
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Appendix B: On the new parameters Fx, CJ\ and c\ 



The Lagrangian (12.11) contains a new field X and three new parameters, namely Fx, o~>i 
and ci, with respect to the usual Lagrangian of Witten, Di Vecchia, Veneziano et al. It is 
therefore natural to ask if the model can be further simplified by simply eliminating some 
parameter. As we have already said, in this paper we are assuming that the parameter 
Fx, which is essentially proportional to the new U(l) axial condensate, is different from 
zero: in section 3 we discuss the relevance of this parameter Fx in the phenomeno logical 
analysis of the strong decays of pseudoscalar mesons. 

Concerning the parameter ui, we cannot say too much. We remind that the usual La- 
grangian of Witten, Di Vecchia, Veneziano et al. is obtained by choosing u>i = 1 (together 
with Fx = 0, i.e., X = 0). At low temperatures one expects that the deviations from 
this Lagrangian are small, in some sense, and therefore one expects that W\ is not much 
different from 1 at low temperatures. On the other side, as already observed in Ref. [6], 
Ui must necessarily be zero when T > T c h, in order to avoid a singular behaviour of the 
anomalous term above the chiral transition: this implies a non trivial behaviour of uj\ with 
the temperature. However, in this paper no particular choice for the value of 0J\(T = 0) 
will be done: it will be considered as a free parameter (apart from the above-mentioned 
limitation for T > T c h). 

The case of the parameter c\ is much more interesting. By putting c\ = 0, i.e., c = [see 
Eq. ill ], into Eqs. (l2~T2D . these reduce to: 



which, when inserted into Eq. (12. lip , lead to the following values for the squared masses 
of the two singlets Si and S2 in the chiral limit: 



Z L 



2A[F- 



- + LF x u\] 



Ql = 0, 



(B.l) 




(B.2) 



The corresponding eigenvectors, written in terms of S n and Sx, are: 





(B.3) 
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Let us observe that Eqs. (IB.3I) and (IB. 21) cannot be derived by simply putting c = into 



Eqs. (I2.14p and (I2.15p . derived in subsection 2.1. This is due to the fact that Eqs. (I2.14p 
and (I2.15P were derived not only assuming that C\ ^ 0, but also taking the large- N limit, 

in which the quantity c = ^= (^t) ( v§) * s ^ ar 9 e ' being of order 0(N) [see Eq. (I2.13P . 
In that case, therefore, one obtains Zl = 0(1) and Ql = 0(1/ N), so that, from Eq. 
(12~TT|) . m\ ~ §^ ~ f3+lf m = ®( l / N ) and & can be identified with the particle r]'. 

Instead, in the particular case in which c\ = [i.e., c = 0], one has that Zl = 0(1/N) 
and Ql = 0, so that, from Eq. (12. lip . Si is massless (in the chiral limit) and therefore it 
does not verify the Witten-Veneziano formula required for the rf '. It is easy to convince 
oneself that, in this particular case c\ = 0, S 2 , having a squared mass of order 0(1/N) 
in the large-iV limit, is just the field which must be identified with the particle rj', as 
required by the Witten-Veneziano mechanism for the solution of the U(l) problem. In 



fact, by virtue of Eqs. (1B.3I) . we can re-write the U(l) axial current jjf), given by Eq. 



(I2.16p . in terms of the fields Si and S 2 : 



5,/i 

where: 

F s 



j£ = -V2Ld, (F Sl Si + F S2 S 2 ) , (B.4) 



F*{ui - 1) + LF 2 x m 
V 'F?(l - uJif + LF^i 



Fs, = , ^ X (B.5) 
y/F^l-coiy + LF^r 

are nothing but the decay constants of the singlet pseudoscalar mesons Si and S 2 , defined 
as: 

(01^(0)1^(^0) = iV2LF Sl p^, 

(0\J^(0)\S 2 (p 2 )) = tV2LF S2 p 2fl . (B.6) 



From Eqs. flB.2|) and (IB.5P one immediately verifies that the field S 2 satifies the Witten- 



Veneziano formula, i.e., 

m s 2 = -#r, ( B - 7 ) 

and, therefore, it is nothing but the field associated with the particle 7/, with a squared 
(non-chiral) mass generated by the anomaly and of order 0(1/N) in the large- iV limit, 
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as required by the Witten-Veneziano mechanism. (Instead, concerning the state Si, even 
if, according to Eqs. flB.6j) and flB.5j) . it is coupled to the U(l) axial current, it is not 



coupled to the topological charge density, i.e., (0\Q(0)\Si(pi)) = ^j-^S'i 777 '! 2 = ^' s i nce ^ 
is massless: therefore it does not appear as an intermediate mesonic state in the spectral 
decomposition of the full topological susceptibility . . .) 

It is interesting to observe that, in this case (differently from the case discussed in subsec- 
tion 2.1), the parameter u>\ plays a fundamental role. In fact, when c\ = 0, the anomalous 
Lagrangian term containing u)\ is the only one which generates a coupling between U and 
X (i.e., between the usual quark-ant iquark pseudoscalar mesons and the exotic singlet 
state). By changing u\ one can "move" the anomaly from U to X. In particular, in the 
case u>i = 1 the anomalous term only depends on U and the field X is decoupled. In 
this case the Lagrangian simply reduces to the sum of the usual Lagrangian written by 
Witten, Di Vecchia, Veneziano et al. for the field U (including the anomalous term) plus 
a non-anomalous Lagrangian for the field X: in this limit the state 5*2, i.e, the rj', reduces 
to the usual quark-antiquark singlet state S n , while the massless state Si reduces to the 
exotic state Sx- On the contrary, in the opposite case u\ = the anomalous term only 
depends on the exotic field X and so the state S2, i.e, the rj', reduces to the exotic state 
Sx, while the massless state Si reduces to the usual quark-antiquark singlet state S n . 
In conclusion, we have found that, in the case in which c\ = 0, in addition to the usual 
L 2 — 1 non-singlet (pseudo-)Goldstone bosons and to the massive singlet S2 = rf, there is 
another singlet Si, which is massless in the chiral limit. This particle is therefore another 
(pseudo-)Goldstone boson which, when including the quark masses, should have a mass 
comparable with that of the other L 2 — 1 non-singlet pseudoscalar mesons. 
In the realistic case L = 3, by diagonalizing the squared mass matrix (12.251) with c = 0, 
we derive the following expressions for the squared masses of Si and S2, at the first order 
in the quark masses and in 1/N: 

2 _ F^l-tm) 2 2 
™ 51 " FZ(l-uir + 3F x ul m °> 
2 /6A\ 2 [2A\ . x2 ^F 2 u 2 2 ,„ \ 

Remembering the definition of m\ = |5(2m + m s ) = |S(m u + m^ + m s ) and Eqs. (I2.23p . 
we immediately see that the squared mass of the singlet Si satifies the following relation: 

m\ = F2n F ' ( |~^L X < < = \{ml + m 2 K± + m 2 KOjRO ) ~ (412 MeV) 2 . (B.9) 
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Even assuming, as already said, that we can identify the singlet S2 with the observed 
singlet rf , no other singlet pseudoscalar meson is observed whose mass satisfies the limit 
(IB. 91) . Our assumption c\ = (together with F x 7^ 0) has thus led us to another U U(1) 
problem". Even if we let C\ be different from zero, but arbitrarily small, i.e., c\ — > 
with all other quantities fixed, since, by virtue of Eqs. fl2.lip - fl2.12l) . the squared masses 
m| 52 in the chiral limit are continuous functions of the parameter ci, we find that 
m 2 Si ~ Ql/Zl — F 2 (i-ui)^ C +LF' i u 2 = ^( Cl ) wn ^ ^ e arbitrarily small and, when including 
quark masses, it will have an upper limit arbitrarily close (from above) to that reported 
in Eq. flB~9l) . 

Therefore, we are forced to discard this possibility (as it leads to wrong predictions for 
the pseudoscalar-meson mass spectrum) and, in the rest of this paper, we shall always 
consider the model in which c\ is different from zero and not too small, so that c = O(N) 
is large. In this case, as we have seen in subsections 2.1 and 2.2, the squared masses of the 
singlet mesons S\ and £2 are given by Eq. ( I2.15P in the chiral limit and by Eqs. (I2.30P 
and (I2.3ip in the realistic case with L = 3 light quark flavours. Therefore, as already said, 
the state S\ has a topological (non-chiral) squared mass of order 0(1/ N) in the large- N 
limit and it is nothing but the particle rj'. Instead, the state S2 is identified with an exotic 
singlet particle rjx, having a large (non-chiral) mass term of order 0(1) in the large-iV 
limit, generated by the (non-zero) coupling constant c\. 
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